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ABSTRACT  The  completed  effort  integrated  the  rigorous  microtevel  (fiber,  matrix, 
sizing)  stress  bound  recently  obtained  by  Lipton  [1.2|  under  AFOSR  sponsorship  with  the 
AFRL-developed  ply  level  multibasis  spline  approximation  stress  analysis  tools  [3-5],  A 
robust  multiscale  analysis  framework  was  developed  and  applied  within  and  beyond  the 
scope  of  the  present  grant.  The  developments  included  extension  of  the  stress  bound 
estimates  in  heterogeneous  materials  with  periodic  microstructure  to  include  residua! 
stress  effects  and,  notably,  regions  of  nonperiodic  microstructure.  Such  regions  included 
ply  interfaces  in  fiber  reinforced  composite  materials,  which  are  present  in  all  laminated 
composite  structures.  The  improved  accuracy  of  the  multiscale  analysis  was  verified  by 
comparison  with  three-dimensional  analysis  of  model  cross  ply  composite  laminate, 
where  fibers  were  modeled  explicitly.  Our  analysis  capability  for  the  first  time  addressed 
the  composition  of  the  repair  patch  of  a  composite  scarf  repair  in  addition  to  more 
conventional  geometric  parameters,  such  as  scarf  angle  and  repair  depth.  We  found  that 
altering  the  conventional  ply  by  ply  replacement  schema  of  the  repair  patch  does  not  lead 
to  premature  failure,  but  may  delay  the  failure  of  the  adhesive  and  increase  the  strength  of 
the  repaired  composite.  This  work  is  well  coordinated  with  the  Materials  Integrity  Branch 
of  the  Air  Force  Research  Laboratory’s  Materials  and  Manufacturing  Directorate,  which 
is  presently  executing  an  experimental  program  to  verify  our  findings.  We  have 
transitioned  the  periodic  RVE  based  micromechanical  failure  criterion  evaluation 
algorithm  to  Sikorsky  Aircraft,  where  it  was  applied  to  failure  analysis  of  helicopter  flex 
beam.  We  are  also  using  our  methodology  to  investigate  the  applicability  of  various 
micromechanical  failure  criteria  to  failure  prediction  of  composites  under  complex  states 
of  stress 
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EXECUTIVE  SUMMARY 

TASK  1.  RVE  models  for  fiber  reinforced  composite  laminates  were  developed 
and  verified  to  provide  accurate  prediction  of  effective  mechanical  and  thermal  expansion 
coefficients  as  a  function  of  fiber  volume  fraction  and  constitutive  properties  of  fiber  and 
matrix.  Anisotropy  of  the  carbon  fiber  was  taken  into  account.  The  quintessential  part  of 
the  task  was  to  demonstrate  the  accuracy  of  application  of  the  simply  coupled  scales 
approximation  (SCSA)  at  the  locations  where  the  microstructure  periodicity  conditions 
were  not  satisfied  such  as  near  free  edges  as  well  as  interfaces  between  plies  with 
different  fiber  orientation.  To  provide  such  estimates,  fully  three-dimensional  models  of 
cross  ply  laminate  section  with  explicit  modeling  of  fiber  and  matrix  were  built.  Figure  I 
illustrates  these  3-D  models.  The  directions  of  the  coordinate  axis  with  the  origin  in  the 
lower  left  corner  and  the  z  -axis  running  in  the  vertical  direction  are  shown  in  Figure  1  as 
well.  The  width  of  the  model  in  the  y-direction  is  equal  to  one  unit  cell  size  In.  The 
boundary  conditions  were  applied  as  follows.  At  the  x=constant  facial  surfaces 
periodicity  conditions  were  applied  along  with  1%  tensile  strain  loading  in  the  x-direction. 
A  displacement  condition  u^O  was  applied  at  the  surface  z=0  and  uy=0  at  the  surface 
y=10/«  so  that  the  3-D  model  corresponds  to  a  quarter  model  of  an  infinitely  long  (in  the 
x-direction)  3  ply  laminate. 

(a)  (b) 
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Figure  I.  3-D  micromechanicat  models  of  (a)  unidirectional/matrix  laminate  and  (b) 
cross-ply  laminates. 
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FIGURE  2.  3-D  model,  SCSA  solution  and  stress  bound  comparison  on  fiber  row 
interfaces 
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The  number  of  unit  cells  modeled  was  determined  by  trial  and  error  by  increasing 
it  until  we  had  enough  cells  in  the  ply  to  achieve  a  good  agreement  between  the  micro 
level  stresses  predicted  by  SCSA  and  the  full  3-D  model  in  the  internal  cells.  We  focused 
on  evaluation  of  the  first  and  second  invariants  of  the  strain  tensor,  which  we  plan  to 
apply  for  failure  prediction  as  proposed  in  [6].  The  results  shown  in  Figure  2  pertain  to 
the  case  of  unidirectional/matrix  ply  laminate  (Figure  I  .a).  Distribution  of  the  first 
invariant  of  the  strain  tensor  as  a  function  of  the  distance  from  the  free  edge  (y- 
coordinate)  is  shown  at  three  different  horizontal  surfaces  at  z=0,  2 In,  3 In.  in  the  middle  of 
the  top  ply  and  away  from  the  free  edge.  Figure  2.a,  the  results  obtained  by  using  the  3-D 
model  and  SCSA  coincide.  Apparently,  the  point  values  of  the  strain  invariants  obtained 
by  SCSA  are  very  accurate  already  for  the  third  unit  cell  away  from  the  free  edge.  The 
bound,  which  appears  to  over  estimate  the  micro  stress  actually  accounts  for  the  strain 
concentrations  on  the  fiber  perimeter,  which  can  also  be  seen  on  the  3-D  solution  if 
plotted  around  the  fibers.  Overall,  the  multiscale  analysis  results  at  the  ply  interface 
significantly  overestimate  the  actual  values  of  the  micro  field  parameters. 

The  parametric  studies  in  [0/90]  s  laminates  revealed  a  less  satisfying  situation. 


FIGURE  3.  Free  edge  problem  statement  (a).  First  strain  invariant  as  a  function  of 
distance  from  the  hole  edge  2.5  fiber  rows  above  interface,  and  (b)  the  second  strain 
invariant  at  the  interface. 

The  distribution  of  the  first  and  second  strain  invariants  in  the  [0/90]s  laminate  loaded  in 
tension  in  the  0(,-direction  is  shown  in  Figure  3.  The  distributions  are  plotted  for  the  fully 
3D  explicit  solution  modeling  each  fiber,  for  the  multi  scale  simply-coupled  scales 
approximation  (SCSA)  solution  and  the  bound.  It  can  be  seen  on  Figure  3,a  that  the 
explicit  and  multiscale  solution  are  in  good  agreement,  even  as  close  as  2.5  fiber 
diameters  from  the  interface.  However,  at  the  interface,  Figure  3.b,  the  multiscale 
solution  and  the  bound  greatly  underestimate  the  amplitude  of  micro-strain  distribution 
resulting  from  exact  solution. 

A  key  achievement  of  the  work  is  the  development  of  the  rigorous 
homogenization  framework  based  on  special  interface  RVE,  shown  in  Figure  4.a.  The 
results  obtained  by  using  the  multiscale 
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FIGURE  4.  Interface  RVE  (a).  The  second  strain  invariant  at  the  interface  as  a  function  of 
the  distance  from  the  free  edge(b). 

solution  based  on  periodic  (dashed  line)  and  interface  RVE  (dotted  line)  are  compared 
with  the  explicit  3D  60  fiber  model  (Figure  lb)»  and  are  displayed  on  Figure  4b.  As  it  can 
be  seen  the  periodic  RVE  based  solution  significantly  underpredicts  the  upper  bound  of  J2 
(second  invariant  of  the  strain  tensor  in  the  matrix  phase)  at  the  interface,  whereas  the 
estimate  provided  by  the  interface  RVE  based  solution  is  much  more  realistic.  These 
results  are  described  in  the  manuscript,  which  we  submitted  to  SIAM  Journal  on 
Multiscale  Analysis. 

TASK  2. 

2.1  Free  edge  dclamination.  A  fully  3-D  ply  level  model  of  special  dog  bone  shaped 
laminated  composite  specimen,  shown  in  Figure  5,  was  developed  according 
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FIGURE  5.  Schematics  of  the  instrumented  dog  bone  type  specimen  for  delamination  onset 
study. 

to  specifications  [7]  to  prediction  delamination  initiation  loads  under  axial  tension 
loading.  Two  levels  of  analysis  are  considered  for  failure  prediction,  ply  stress  based 
failure  criterion  and  micro-level  failure  criterion. 

A  comparison  of  ply  level  stress  distributions  in  laminates  with  various 
layups  is  shown  on  Figure  6.  The  ply  orientations  in  the  laminates  were  [0/-0/0/-0/9O]s 
for  0=15°,  201’,  25l1,  30*’,  3511  and  40(l,  The  coordinate  direction  and  origin  are  shown  in 
Figure  5.  Tensile  loading  is  applied  in  the  x-direction.  According  to  the  experimental 
study,  the  delamination  initiation  load  for  the  laminates  decreased  by  increasing  0.  We 
examined  the  x=0  cross  section  of  the  laminate,  which  is  in  the  gage  section,  where  the 
delamination  onset  is  monitored.  Figure  6  shows  the  distribution  of  transverse  normal  and 
shear  stress  components  and  cryz  at  the  surfaces  where  they  exhibited  the  maximum 
values,  which  w'ere  the  midsurface  and  0/-0  surface  for  two  laminates  with  0=15°  and 
0=40°. 


a) 


b) 
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FIGURE  6.  Stress  component  comparison  of  the  ply  level  solution  between  laminates 
with  0=40°  (open  markers)  and  0=15°  laminates  (filled  markers):  (a)  normal  cr^ stress  at 
the  mid-surface,  z=0  and  (b)  transverse  shear  stress  <ryz  at  the  0/-0  interface  (b). 

Both  laminates  are  loaded  with  the  same  average  strain.  As  seen  in  Figure  6  it  is  the 
transverse  normal  stress  <ja  which  dominated  over  the  shear  stress  magnitude  in  both 
laminates  and  its  level  in  the  0=4O()  laminate  is  higher  near  the  free  edge  then  in  the  0=15° 
laminate.  This  correlated  well  with  the  experimental  data  showing  higher  delamination 
initiation  loads  in  the  0=1 511  plies. 

Multiscale  analysis  solution  and  calculated  bounds  for  the  first  and  second 
invariants  of  the  strain  tensor  in  the  matrix  phase  are  shown  on  Figure  7.  In  the  example 
at  hand  we  are  less  interested  in  fiber  characteristics,  since  the  initial  failure  mode  is 
matrix  cracking  and  delamination.  Figure  7  displays  the  values  of  the  first  and  second 
invariants  of  the  strain  tensor  calculated  at  the  midsurface  in  the  x=0  cross  section  as  a 
function  of  distance  from  the  hole  edge.  Figure  7.a  displays  results  obtained  by  using  ply 
level  analysis  and  Figure  7.b  results  obtained  by  using  SCSA  for  the  two  laminates.  In 
this  case  both  solutions  yield  qualitatively  similar  results.  Significantly  larger  values  of 
the  micro  level  stress  bounds  reflect  the  amplitude  of  the  strain  oscillation  in  the  unit  cell. 
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Figure  7.  Comparison  of  the  distributions  of  strain  tensor  invariants  predicted  in  the 
matrix  phase  at  the  midsurface  at  x=0  as  a  function  of  distance  from  the  free  edge  by 
using  the  pty  level  solution  (a)  and  (b)  the  bounds  of  strain  tensor  invariants  predicted  by 
SCSA. 

2,2.  Normal  open-hole.  In  this  task,  we  compared  the  distribution  of  the  level  sets  of 
matrix  phase  strain  invariants  calculated  by  using  the  multiscale  solution  and  the  level 
sets  of  the  strain  invariants  calculated  from  the  macro  level  strains.  We  also  compared  the 
normalized  matrix  phase  strain  invariants  calculated  by  using  the  multiscale  solution  and 
the  level  sets  of  normalized  macro-level  stresses  and  their  combinations.  Such  a 
comparison  is  motivated  by  the  existing  approach  to  failure  prediction  in  composite 
laminates.  Thus,  the  failure  criteria  applied  on  the  micro  level  represent  the  Strain 
Invariant  Failure  Theory  (SIFT)  [6]  and  consist  of  the  first  and  second  invariants  of  the 
strain  tensor.  On  the  ply  level,  however,  the  most  common  failure  criteria  are  applied  in 
the  stress  space.  We  have  calibrated  both  failure  criteria  with  the  same  experimental 
values. 
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Comparison  of  invariants  of  the  homogenized  strain  field  to  those  of  the  macro  level 
strain  invariants.  The  first  invariant  bound  and  second  invariant  bound  will  be  denoted 
£m((?(i/1))  and  nM(e(«°)),  respectively.  A  0°  unidirectional,  4  in.  x  1  in.  rectangular  plate 
with  a  0.25  in.  diameter  open  circular  hole  was  considered.  Ply  thickness  was  0.005  in. 
A  1%  tension  strain  in  the  ^-direction  was  applied  as  a  displacement  condition  («x  =  0,04) 
on  the  surface  x  =  2.  In  addition,  the  ^-displacement  was  restricted  on  the  surface  x  =  2. 
Similarly,  the  ^-displacement  and  y-displacement  were  restricted  on  the  surface  *  =  -2, 
Midplane  symmetry  was  provided  for  the  laminate  by  the  boundary  condition  u:  =  0  on 
the  surface  z  =  0,  Thus  the  results  correspond  to  a  2-ply  model. 

Figure  8(a)  and  8(b)  display  the  level  sets  of  the  second  invariant  of  the 
homogenized  strain  field  and  the  bounds  from  the  multiscale  solution  nM(e(«0)), 
respectively.  Figure  8  also  suggests  that  the  actual  strain  within  the  composite  can  be 
much  larger  than  the  state  predicted  by  the  homogenized  strain  field. 


(a)  Invariant  of  the  homogenized  strain  field  (b)  Macroscopic  failure  criteria  nM(e(w0)) 


Figure  8:  Second  strain  invariant  comparison  for  unidirectional  0°  laminate 

Effect  of  thermal  prestress.  The  model  for  the  thermal  prestress  problem  was  identical 
to  the  model  of  the  previous  problem,  except  that  we  used  a  4-ply  laminate  (the  z- 
symmetry  condition  corresponds  to  an  8-ply  symmetric  laminate).  Each  ply  in  the 
laminate  had  a  thickness  of  0.005  in.,  and  thus  the  8-ply  specimen  was  0.04  in.  thick  in 
the  z-direction.  The  8-ply  laminate  was  an  IM7/977-3  [0/45/-45/90]s  composite.  The 
symmetry  plane  was  the  surface  z  =  0.  The  strain  invariant  bounds  in  the  matrix  phase 
were  computed  at  four  z-heights  in  the  model,  corresponding  to  the  middle  (with  respect 
to  2)  of  each  ply.  The  first  strain  invariant  in  the  45°  ply  is  shown  in  Figure  9.  The 
prestress  appears  to  increase  the  invariant  value  by  one  order  of  magnitude  on  nearly  the 
entire  ply.  This  trend  is  representative  of  the  pre stress  effects  in  the  0°  and  ±45°  plies. 
The  pre  stress  did  increase  the  first  invariant  value  for  the  90°  ply,  but  by  only  a  few 
tenths. 

The  values  of  the  bounds  for  the  second  nM(e(«0))  in  the  45°  ply  are  displayed  in 
Figure  10.  It  is  again  clear  from  the  figures  that  the  prestress  increases  the  second 
invariant  values  on  nearly  the  entire  plate.  The  prestress  does  not,  however,  appear  to 
significantly  increase  the  maximum  of  the  second  invariant  inside  the  ply.  It  is  seen  that 
the  prestress  increases  the  volume  of  the  sets  experiencing  the  higher  strain  states.  This 
volume  can  be  crucial  in  failure  predictions,  and  thus  its  increase  is  significant.  It  is 
important  to  note  that  second  invariant  values  in  the  ±45°  plies  are  about  twice  as  large  as 
the  values  in  the  0°  and  90°  plies.  The  off-axis  fiber  orientation  angles  clearly  experience 
larger  shear  strains. 
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It  is  important  to  note  that  the  largest  strains  in  each  ply  are  concentrated  very 
near  the  hole  edge.  Hence  contours  approaching  these  extreme  values  are  not  visible  in 
these  figures.  The  invariant  contours  plotted  show  the  general  trends  of  the  invariants  in 
each  ply,  but  do  not  represent  the  extreme  invariant  values  that  occur.  For  the  first 
invariant,  the  largest  level-line  plotted  in  each  figure  is  roughly  50%  of  the  maximum 
invariant  in  the  ply.  For  the  second  invariant,  the  largest  level-line  shown  in  each  set  of 
figures  is  about  20%  of  the  maximum  invariant  in  the  ply. 


(a)  jCm(£>(m0))  without  prestress  (b)  XM(e(w0))  with  prestress 


Figure  9:  Macroscopic  failure  criteria  £M(e(u0))  in  the  +45°  ply 


(a)  nM(e(w0))  without  prestress  (b)  nM(e(u0))  with  prestress 


Figure  10:  Macroscopic  failure  criteria  nM(e(u0))  in  the  +45°  ply 


Comparison  of  cubic  RVE  and  hexagonal  RVE.  We  created  a  hexagonal  unit  cell 
model.  Figure  5  shows  a  schematic  of  the  cubic  and  hexagonal  unit  cell  models, 
respectively.  The  macroscopic  computer  model  for  this  problem  was  an  8-ply,  symmetric 
laminate  with  an  open  hole,  and  was  identical  in  dimension  and  loading  to  the  model  for 
the  previous  problem.  Data  was  taken  with  and  without  prestress,  and  the  bounds  of  the 


(a)  Cubic  Unit  Cell  (b)  Hexagonal  Unit  Cell 
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FIGURE  1 1 :  Schematic  RVE  models  of  two  fiber  reinforced  geometries 

Figure  12  shows  the  first  invariant  bounds  £M(e(a°))  with  and  without  prestress. 
The  Figure  suggests  once  again  that  the  prestress  increases  the  maximum  value  of  the 
macroscopic  failure  criteria,  as  well  as  increasing  its  value  on  the  entire  plate. 
Comparing  Figure  9  and  Figure  12,  it  is  also  evident  that  the  strain  state  is  smaller  when 
the  microstructurc  is  hexagonal ly  periodic. 


(a)  £M(e(w0))  without  prestress  (b)  £M(e(w0))  with  prestress 


FIGURE  1 2:  Bounds  for  the  first  strain  invariant  £M(e(w11))  in  the  +45°  ply 


We  note  here  that  the  color  scales  are  the  same  for  Figure  9  and  Figure  12,  as  well 
as  for  Figure  10  and  Figure  13.  Thus  the  corresponding  figures  for  the  cubic  RVE  and 
the  hexagonal  RVE  can  be  compared  directly  by  observing  the  color.  Figure  13  plots 
second  strain  invariant  bound  nM(e(w0))  with  and  without  prestress.  Once  again,  the 
macroscopic  failure  criteria  values  calculated  using  the  hexagonal  RVE  were  less  than  the 
values  calculated  using  the  cubic  RVE. 


(a)  nM(e(«0))  without  prestress 


Figure  13:  The  bound  for  the  second  strain  invariant  nM(tJ(w0))  in  the  +45°  ply 


Failure  analysis  of  a  symmetric  laminate  with  an  open  hole.  We  compared  the 
normalized  bounds  of  the  matrix  phase  strain  invariants  with  the  normalized  macro-level 
stress  values  and/or  their  combinations.  Experimental  strength  data  were  used  to  back- 
calculale  critical  strain  invariant  values  for  the  matrix  phase  of  the  composite.  These 
critical  values  were  used  to  normalize  the  strain  invariant  bounds.  Hashin  criteria  and 
maximum  stress  criteria  were  applied  to  the  ply  level  (homogenized)  stresses.  The  over¬ 
stressed  (over-strained)  zone  is  the  volume  of  the  composite  in  which  the  value  of  the 
failure  criteria  meets  or  exceeds  any  of  the  critical  strength  parameters  (indicated  by 
yellow  in  Figures  15  -  19). 

The  computer  model  for  this  experiment  was  identical  in  dimension  to  the  model 
of  the  previous  two  experiments.  The  loading,  however,  was  reduced.  Since  complete 
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loss  of  load  carrying  capacity  typically  occurs  between  0.6%  strain  and  0.7%  strain,  we 
computed  over-stressed  zones  at  0.3%  and  0.6%  uniaxial  tension  (applied  via 
displacement  on  the  surface  x  =  2).  All  other  boundary  conditions  were  the  same  as  in 
the  prior  experiments. 

Figures  14,  15,  and  16  display  the  overstressed  zones  for  the  phase  strain  invariant, 
Hashin  failure  criteria,  and  maximum  stress  failure  criteria,  respectively  at  0.3%  tension 
strain  loading.  At  such  a  small  loading,  the  0°  ply  has  very  small  overstressed  zones  just 
at  the  hole  surface.  The  45°  ply  shows  moderately  smalt  over-stressed  zones  protruding 
from  the  hole  surface  out  toward  the  ply  edge.  The  -45°  ply  (not  shown)  was  similar  to 
the  45°  ply,  and  the  90°  ply  (also  not  shown)  was  almost  completely  overstressed  at  just 
0.3%  strain.  The  three  failure  theories  agree  quite  well  on  the  location  and  size  of  the 
over-stressed  zones.  Moreover,  the  theories  appear  to  agree  on  the  failure  mode.  The 
dashed-line  boundary  around  the  macroscopic  failure  criteria  overstrained  zone  (Figure 
14(a))  indicates  shear  failure.  The  solid-line  boundary  around  the  Hashin  overstressed 
zone  (Figure  13(a))  indicates  failure  via  a  quadratic  combination  of  transverse  tension 
and  in-plane  shear  forces.  The  boundary  around  the  maximum  stress  over-stressed  zone 
(Figure  14(a))  indicates  failure  by  the  in-plane  shear  force.  Thus  the  results  suggest  that 
the  three  theories  agree  that  the  shear  forces  will  initiate  failure  at  the  hole  edge  in  the  0° 

ply- 

Figures  17,  18,  and  19  visualize  the  over-stressed  zones  for  the  micro  level  strain 
bounds,  Hashin  failure  criteria,  and  max  stress  failure  criteria,  respectively  at  0.6% 
tension  strain  loading.  The  90°  and  ±45°  plies  all  have  large  over-stressed  zones  (only 
+45°  ply  is  shown).  The  small  overstressed  zones  very  near  the  hole  edge  in  the  0°  ply 
are  the  shear  overstressed  zones  seen  at  0.3%  strain.  These  shear  stresses  near  the  hole 
edge  produce  the  characteristic  cracks  running  parallel  to  the  fibers  in  the  0°  ply.  The 
slightly  larger  overstressed  zones  in  the  0°  ply  appear  to  be  transverse  tension  failure. 
These  overstressed  zones  are  believed  to  be  an  effect  of  the  ply  stacking  sequence.  Once 
again,  the  micro-level  strain  invariants  based  failure  criteria  appear  to  be  consistent  with 
the  predictions  of  the  Hashin  failure  criteria  and  the  maximum  stress  failure  criteria.  In 
addition  to  the  three  different  failure  criteria  being  in  agreement  on  the  size  and  location 
of  the  overstressed  zones,  they  also  once  again  appear  to  agree  on  the  failure  mode  in 
each  case. 
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(a)  0°  Ply 


(b)  +45°  Ply 


FIGURE  14:  Micro  level  second  strain  invariant  overstrained  zone  at  0.3%  strain 


(a)  0°  Ply 


(b)  +45°  Ply 


FIGURE  15:  !  lashin  failure  criteria  overstressed  zone  at  0.3%  strain 


(a)  0°  Ply 


(b)  +45°  Ply 


Figure  16:  Maximum  stress  failure  criteria  overstressed  zone  at  0.3%  strain 


(a)  0°  Ply 


(b)  +45°  Ply 


Figure  17:  Micro-level  second  invariant  overstrained  zone  at  0.6%  strain 
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(a)  0°  Ply  (b)  +45°  Ply 

Figure  18:  Hashin  failure  criteria  overstressed  zone  at  0.6%  strain 


(a)  0°  Ply 


(b)  +45°  Ply 


Figure  19:  Maximum  stress  failure  criteria  overstressed  zone  at  0.6%  strain 


TASK  3.  In  this  task  we  tried  to  understand  the  failure  sequence  of  the  composite  scarf 
repair  patch  and  construct  a  meaningful  goal  function  for  strength  retention  optimization. 
The  main  development  in  this  task  was  to  model  nonlinear  deformation  and  failure  in  the 
adhesive.  The  deformation  of  the  composite  and  the  patch  was  modeled  with  less  detail 
on  ply  level.  We  heavily  utilized  experimental  data  obtained  by  Capt.  B.M.  Cook  [8], 
which  included  fabrication  and  testing  of  the  scarf  repair  as  shown  on  Figure  20. 


FIGURE  20.  Experimental  setup  for  tensile  testing  (a)  and  a  plane  view  of  the  failed  scarf 
repaired  specimen  (b). 

The  key  finding  in  this  task  was  that  modeling  of  the  failure  of  the  adhesive  in  the 
nonlinear  regime  is  crucial  to  predicting  the  strength  of  the  repair.  We  found  that  despite 
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the  low  stiffness  the  adhesive  transfers  enough  load  in  the  elastic  regime  to  practically 
remove  the  stress  concentration  due  to  scarfing  and  patching  thus  leading  to  95%  strength 
retention  prediction.  It  is  the  adhesive  failure  that  drastically  increases  the  stress 
concentration  in  the  parent  panel  and  leads  to  premature  failure.  Spline  approximation 
based  3D  ply  level  analysis  was  used  for  failure  prediction.  The  plate,  adhesive  and  the 
repair  patch  models  are  shown  on  Figure  21. 

Patch 


Figure  2 1 .  Computational  model  for  the  scarfed  plate,  patch  and  adhesive. 

A  nonlinear  elastic  model  was  utilized  to  model  the  adhesive  At  the  same  time  the  fiber 
failure  in  the  scarfed  panel  was  predicted  by  using  the  critical  failure  volume  (CFV) 
method  recently  proposed  in  [9-10].  The  average  applied  load  for  fiber  failure  was 
predicted  at  each  load  level  characterized  by  a  state  of  delamination  in  the 
adhesive/scarfed  panel  interface.  Initially,  the  stress  distribution  in  the  repaired  panel  is 
uniform  and  the  fiber  failure  stress  predicted  is  identical  to  that  of  virgin  laminate.  The 
increase  of  the  applied  load  leads  to  initiation  and  growth  of  the  delamination  between 
the  repair  patch  and  the  plate.  The  stress  distribution  in  the  patch  becomes  non  uni  form 
with  the  stress  concentration  typical  of  open  hole  composites.  At  some  applied  load  the 
average  predicted  fiber  failure  load  and  the  applied  load  become  equal.  This  load  value  is 
considered  the  failure  load  for  the  repaired  panel. 

The  unidirectional  properties  of  the  IM6/3501-6  material  were  as  follows: 
En=175.27Gpa,  Ej2=E33=9.79Gpa,  v 1 3=  V| 3=0.33,  V23=0.49,  Gi2=G  1 3=5.5  IGpa  and 
G23=2.96Gpa,  with  the  coefficient  of  thermal  expansion  of  an=4.!0'6C°  and  ct22“32  10' 
6C°  The  Weibull  parameters  for  the  strength  in  the  fiber  direction  were  equal  to 
Xi=2.06Gpa.  Vo=l  638.7mm3  and  a=40.  The  value  of  the  Weibull  modulus  was  taken 
from  Wisnom  et  al.  [11],  The  /mjn=0.266mm  was  calculated  according  to  reference  [9-10]. 
The  initial  elastic  properties  of  the  adhesive  were  E=3.1Gpa,  v=0.38  and  the  thermal 
expansion  coefficient  was  a=0.7  10‘6C°.  The  ply  and  adhesive  thickness  were  h=0.13mm. 
The  nonlinear  shear  deformation  curve  i(y)  for  the  FM-300M  adhesive  was  provided  by 
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Structural  Integrity  Branch  of  the  Air  Force  Research  Laboratory’s  Materials  and 
Manufacturing  Directorate  [12]  and  is  shown  in  Figure  21.  Three-  dimensional  nonlinear 
deformation  of  the  adhesive  was  describer  by  assuming  nonlinear  elasticity  with  constant 

poisson  ratio  and  shear  modulus  equal  to  G=r(2/  V3  sVm),  where  zvm~  is  Von-Mises  strain, 
where  G(y)=dx(y)/dy.  The  temperature  change  AT=-155°C. 


Figure  2 1 .  Nonlinear  shear  deformation  curve  of  the  for  the  FM-300M  adhesive. 


At  low  loads  prior  to  the  adhesive  failure  the  stress  Held  in  the  scarfed  panel  is  quite 
uniform,  as  shown  on  the  Figure  22.  No  stress  concentration  typical  of  an  open  hole  is 
present.  The  average  fiber  failure  load  is  predicted  equal  to  that  of  virgin  panels. 

(a)  b) 


Figure  22.  Normal  crxx  stress  distribution  (a)  in  the  adhesive  and  (b)  the  scarfed  plate  prior 
to  any  delamination  on  the  adhesive/panel  interface 
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In  the  case  of  scarfed  open  hole  panels  a  severe  stress  concentration  typical  of  an  open 
hole  exists,  Figure  23. a.  The  strength  of  such  panels  is  about  1/3  of  the  virgin  panels  and 
is  predicted  accurately  by  using  the  CFV  method  [9-10],  Thus  the  strength  of  the  repair  is 
determined  by  the  adhesive  failure  and/or  its  nonlinear  behavior.  The  stress  distribution 
in  the  parent  panel  corresponding  to  partially  failed  adhesive  is  shown  on  Figure  23b. 
This  stress  concentration  is  less  severe  than  in  the  case  of  open  scarfed  hole  and  the 
predicted  strength  values  of  the  repaired  panels  are  at  the  level  of  2/3  of  the  initial  virgin 
panel  strength. 


(a)  b) 


Figure  23.  Normal  cr**  stress  distribution  in  the  scarfed  and  repaired  panels,  (a)  Open 
scarfed  hole  and  (b)  partially  delaminated  adhesive  at  90%  of  the  scarf  failure  load. 


Experimental  data  and  predicted  strength  values  for  small  and  large  size  unnotched 
tensile  coupons,  as  well  as  strength  of  scarfed  open  hole  and  scarfed  repaired  panels  are 
shown  on  Figure  24.  Significant  reduction  in  experimentally  measured  strength  of  large 
virgin  panels  as  compared  to  small  coupons  is  attributed  to  grip  failure  of  most  large 
virgin  specimens.  Only  5%  strength  reduction  due  to  volume  increase  is  predicted 
theoretically.  All  unrepaired  and  repaired  scarfed  laminates  exhibited  brittle  failure 
through  the  cavity,  thus  rendering  the  test  data  valid  for  comparison  with  theoretical 


Figure  24.  Experimental  and  numerical  results 
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experimentally. 

TASK  4.  Development  of  an  automated  mutiscale  mechanism  within  the  B-splinc 
Approximation  Method  (BSAM)  computer  program  was  accomplished  ahead  of  schedule 
and  utilized  to  produce  all  results  demonstrated  in  this  report. 

TASK  5. 

In  the  course  of  the  project,  we  have  understood  that  our  original  vision  of  repair 
optimization  based  on  elastic  properties  of  adhesive  was  too  simplistic.  Indeed  as  seen 
from  the  previous  task  the  repair  strength  is  a  direct  function  of  the  adhesive  failure, 
which  determines  the  load  transfer  between  the  repair  patch  and  the  parent  laminate.  An 
optimization  algorithm  which  would  be  based  on  a  repair  strength  as  a  goal  function  will 
be  highly  impractical,  requiring  nonlinear  analysis  at  every  optimization  iteration  and 
preventing  the  use  of  gradient  methods  [13]  for  optimization.  These  complications 
delayed  our  progress  and  required  to  propose  alternative  optimization  ideas.  The  one  we 
have  selected  for  implementation  is  to  minimize  an  elastic  Lpnorm  of  a  stress  invariant 
failure  criteria  in  the  adhesive.  Such  an  approach  is  targeted  to  the  delay  of  failure 
initiation  in  the  adhesive.  Although  not  completed,  we  have  obtained  encouraging  results. 
Recall  Figure  24,  where  the  predicted  strength  of  the  repaired  composite  was  at  60Ksi. 
This  result  was  obtained  for  a  traditional  ply  by  ply  replacement  composition  of  the  repair 
patch.  After  we  switched  two  ply  orientations,  namely  0"  and  90°  inside  this  patch,  we 
delayed  the  onset  of  adhesive  failure  to  65Ksi.  We  expect  to  further  delay  the  adhesive 
failure  by  application  of  the  optimization  algorithm.  These  developments  are  being 
closely  followed  by  the  Structural  Integrity  Branch  of  the  Air  Force  Research 
Laboratory's  Materials  and  Manufacturing  Directorate,  which  is  presently  conducting  an 
experimental  program  to  verify  our  findings, 

CONCLUSIONS 

1)  We  integrated  the  rigorous  microleve!  (fiber,  matrix)  stress  bound  estimates  by 
using  simply  coupled  scale  analysis,  with  the  spline  approximation  stress  analysis  tools  - 
BSAM.  A  robust  multiscale  analysis  framework  was  developed  and  applied  within  and 
beyond  the  scope  of  the  present  grant. 

2)  A  key  development  was  the  extension  of  the  stress  bound  estimates  in 
heterogeneous  materials  with  periodic  microstructure  to  include  regions  of  nonperiodic 
microstructurc.  Such  regions  included  ply  interfaces  in  fiber  reinforced  composite 
materials,  which  are  present  in  all  laminated  composite  structures.  The  improved 
accuracy  of  the  multiscale  analysis  was  verified  by  comparison  with  three-dimensional 
analysis  of  model  cross  ply  composite  laminate,  where  fibers  were  modeled  explicitly. 

3)  The  developed  capabilities  were  applied  to  composite  scarf  repair  optimization. 
Our  analysis  addressed,  for  the  first  time,  the  composition  of  the  repair  patch  in  addition 
to  more  conventional  geometric  parameters,  such  as  scarf  angle  and  repair  depth.  We 
found  that  altering  the  conventional  ply  by  ply  replacement  schema  of  the  repair  patch 
does  not  lead  to  premature  failure,  but  delays  the  failure  of  the  adhesive  and  increases  the 
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strength  of  the  repaired  composite.  This  work  is  coordinated  with  Structural  Integrity 
Branch  of  the  Air  Force  Research  Laboratory’s  Materials  and  Manufacturing  Directorate, 
which  is  presently  executing  an  experimental  program  to  verify  our  findings. 

4)  We  have  transitioned  the  periodic  RVE  based  micromechanical  failure  criterion 
evaluation  algorithm  to  Sikorsky  Aircraft,  where  it  was  applied  to  failure  analysis  of  a 
helicopter  flex  beam.  We  are  also  applying  our  methodology  to  investigate  applicability 
of  various  micromechanica]  failure  criteria  to  failure  prediction  of  composites  under 
complex  states  of  stress. 
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ABSTRACT  The  completed  effort  integrated  the  rigorous  microlevel  (fiber,  matrix, 
sizing)  stress  bound  recently  obtained  by  Lipton  [1,2]  under  AFOSR  sponsorship  with  the 
AFRL-developed  ply  level  multibasis  spline  approximation  stress  analysis  tools  [3-5],  A 
robust  multiscale  analysis  framework  was  developed  and  applied  within  and  beyond  the 
scope  of  the  present  grant.  The  developments  included  extension  of  the  stress  bound 
estimates  in  heterogeneous  materials  with  periodic  microstructure  to  include  residual 
stress  effects  and,  notably,  regions  of  nonperiodic  microstructure.  Such  regions  included 
ply  interfaces  in  fiber  reinforced  composite  materials,  which  are  present  in  all  laminated 
composite  structures.  The  improved  accuracy  of  the  muiliscaie  analysis  was  verified  by 
comparison  with  three-dimensional  analysis  of  model  cross  ply  composite  laminate, 
where  fibers  were  modeled  explicitly.  Our  analysis  capability  for  the  first  time  addressed 
the  composition  of  the  repair  patch  of  a  composite  scarf  repair  in  addition  to  more 
conventional  geometric  parameters,  such  as  scarf  angle  and  repair  depth.  We  found  that 
altering  the  conventional  ply  by  ply  replacement  schema  of  the  repair  patch  does  not  lead 
to  premature  failure,  but  may  delay  the  failure  of  the  adhesive  and  increase  the  strength  of 
the  repaired  composite.  This  work  is  well  coordinated  with  the  Materials  Integrity  Branch 
of  the  Air  Force  Research  Laboratory's  Materials  and  Manufacturing  Directorate,  which 
is  presently  executing  an  experimental  program  to  verify  our  findings.  We  have 
transitioned  the  periodic  RVE  based  micromechanical  failure  criterion  evaluation 
algorithm  to  Sikorsky  Aircraft,  where  it  was  applied  to  failure  analysis  of  helicopter  flex 
beam.  We  are  also  using  our  methodology  to  investigate  the  applicability  of  various 
micromechanical  failure  criteria  to  failure  prediction  of  composites  under  complex  states 
of  stress 
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EXECUTIVE  SUMMARY 

TASK  1.  RVE  models  for  fiber  reinforced  composite  laminates  were  developed 
and  verified  to  provide  accurate  prediction  of  effective  mechanical  and  thermal  expansion 
coefficients  as  a  function  of  fiber  volume  fraction  and  constitutive  properties  of  fiber  and 
matrix.  Anisotropy  of  the  carbon  fiber  was  taken  into  account.  The  quintessential  part  of 
the  task  was  to  demonstrate  the  accuracy  of  application  of  the  simply  coupled  scales 
approximation  (SCSA)  at  the  locations  where  the  microstructure  periodicity  conditions 
were  not  satisfied  such  as  near  free  edges  as  well  as  interfaces  between  plies  with 
different  fiber  orientation.  To  provide  such  estimates,  fully  three-dimensional  models  of 
cross  ply  laminate  section  with  explicit  modeling  of  fiber  and  matrix  were  built.  Figure  1 
illustrates  these  3-D  models.  The  directions  of  the  coordinate  axis  with  the  origin  in  the 
lower  left  corner  and  the  z  -axis  running  in  the  vertical  direction  are  shown  in  Figure  l  as 
well.  The  width  of  the  model  in  the  y-direction  is  equal  to  one  unit  cell  size  to.  The 
boundary  conditions  were  applied  as  follows.  At  the  x=constant  facial  surfaces 
periodicity  conditions  were  applied  along  with  1%  tensile  strain  loading  in  the  x-direction, 
A  displacement  condition  uz=0  was  applied  at  the  surface  z=0  and  uy=0  at  the  surface 
y=l0/n  so  that  the  3-D  model  corresponds  to  a  quarter  mode!  of  an  infinitely  long  (in  the 
x-direction)  3  ply  laminate. 

(a)  (b) 
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Figure  1. 3-D  micromechanical  models  of  (a)  unidirectional/mairix  laminate  and  (b) 
cross-ply  laminates, 
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Figure  2.  3-D  model,  SCSA  solution  and  stress  bound  comparison  on  fiber  row 
interfaces 
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The  number  of  unil  cells  modeled  was  determined  by  trial  and  error  by  increasing 
it  until  we  had  enough  cells  in  the  ply  to  achieve  a  good  agreement  between  the  micro 
level  stresses  predicted  by  SCSA  and  the  full  3-D  model  in  the  interna!  cells.  We  focused 
on  evaluation  of  the  first  and  second  invariants  of  the  strain  tensor,  which  we  plan  to 
apply  for  failure  prediction  as  proposed  in  [6].  The  results  shown  in  Figure  2  pertain  to 
the  case  of  unidirectional/matrix  ply  laminate  (Figure  l.a).  Distribution  of  the  first 
invariant  of  the  strain  tensor  as  a  function  of  the  distance  from  the  free  edge  (y- 
coordinate)  is  shown  at  three  different  horizontal  surfaces  at  z=0,  2 lo.  3 /«.  In  the  middle  of 
the  top  ply  and  away  from  the  free  edge,  Figure  2. a.  the  results  obtained  by  using  the  3-D 
model  and  SCSA  coincide.  Apparently,  the  point  values  of  the  strain  invariants  obtained 
by  SCSA  are  very  accurate  already  for  the  third  unit  cell  away  from  the  free  edge.  The 
bound,  which  appears  to  over  estimate  the  micro  stress  actually  accounts  for  the  strain 
concentrations  on  the  fiber  perimeter,  which  can  also  be  seen  on  the  3-D  solution  if 
plotted  around  the  fibers.  Overall,  the  multiscale  analysis  results  at  the  ply  interface 
significantly  overestimate  the  actual  values  of  the  micro  field  parameters. 

The  parametric  studies  in  [0/90]s  laminates  revealed  a  less  satisfying  situation. 


FIGURE  3.  Free  edge  problem  statement  (a).  First  strain  invariant  as  a  function  of 
distance  from  the  hole  edge  2.5  fiber  rows  above  interface,  and  (b)  the  second  strain 
invariant  at  the  interface. 

The  distribution  of  the  first  and  second  strain  invariants  in  the  [0/90]s  laminate  loaded  in 
tension  in  the  O(1-direction  is  shown  in  Figure  3.  The  distributions  are  plotted  for  the  fully 
3D  explicit  solution  modeling  each  fiber,  for  the  multi  scale  simply-coupled  scales 
approximation  (SCSA)  solution  and  the  bound,  it  can  be  seen  on  Figure  3. a  that  the 
explicit  and  multiscale  solution  are  in  good  agreement,  even  as  close  as  2.5  fiber 
diameters  from  the  interface.  However,  at  the  interface,  Figure  3.b,  the  multiscale 
solution  and  the  bound  greatly  underestimate  the  amplitude  of  micro-strain  distribution 
resulting  from  exact  solution. 

A  key  achievement  of  the  work  is  the  development  of  the  rigorous 
homogenization  framework  based  on  special  interface  RVE,  shown  in  Figure  4. a.  The 
results  obtained  by  using  the  multiscale 
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FlGURl-  4.  Interface  RVE  (a).  The  second  strain  invariant  at  the  interface  as  a  function  of 
the  distance  from  the  free  edge(b). 

solution  based  on  periodic  (dashed  line)  and  interface  RVE  (dotted  line)  are  compared 
with  the  explicit  3D  60  fiber  model  (Figure  lb),  and  are  displayed  on  Figure  4b.  As  it  can 
be  seen  the  periodic  RVE  based  solution  significantly  underpredicts  the  upper  bound  of  J; 
(second  invariant  of  the  strain  tensor  in  the  matrix  phase)  at  the  interface,  whereas  the 
estimate  provided  by  the  interface  RVE  based  solution  is  much  more  realistic.  These 
results  are  described  in  the  manuscript,  which  we  submitted  to  SIAM  Journal  on 
Multiscale  Analysis. 

TASK  2. 

2.1  Free  edge  delami nation.  A  fully  3-D  ply  level  model  of  special  dog  bone  shaped 
laminated  composite  specimen,  shown  in  Figure  5,  was  developed  according 
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FIGURE  5.  Schematics  of  the  instrumented  dog  bone  type  specimen  for  delamination  onset 
study. 

to  specifications  [7]  to  prediction  delamination  initiation  loads  under  axial  tension 
loading.  Two  levels  of  analysis  are  considered  for  failure  prediction,  ply  stress  based 
failure  criterion  and  micro-level  failure  criterion. 

A  comparison  of  ply  level  stress  distributions  in  laminates  with  various 
layups  is  shown  on  Figure  6.  The  ply  orientations  in  the  laminates  were  [9/-0/0/-0/9O]s 
for  0=15°,  20°,  25°,  30°,  35IJ  and  40°.  The  coordinate  direction  and  origin  are  shown  in 
Figure  5.  Tensile  loading  is  applied  in  the  x-direetion.  According  to  the  experimental 
study,  the  delamination  initiation  load  for  the  laminates  decreased  by  increasing  0.  We 
examined  the  x=0  cross  section  of  the  laminate,  which  is  in  the  gage  section,  where  the 
delamination  onset  is  monitored.  Figure  6  shows  the  distribution  of  transverse  normal  and 
shear  stress  components  a?*  and  cryz  at  the  surfaces  where  they  exhibited  the  maximum 
values,  which  were  the  midsurface  and  0/-0  surface  for  two  laminates  with  9=1 5l>  and 
0=40°. 


a) 


b) 
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FIGURE  6.  Stress  component  comparison  of  the  ply  level  solution  between  laminates 
with  0=40°  (open  markers)  and  0=1 5()  laminates  (filled  markers):  (a)  normal  stress  at 
the  mid-surface,  z=0  and  (b)  transverse  shear  stress  cyi-  at  the  0/-0  interface  (b). 

Both  laminates  are  loaded  with  the  same  average  strain.  As  seen  in  Figure  6  it  is  the 
transverse  normal  stress  which  dominated  over  the  shear  stress  magnitude  in  both 
laminates  and  its  level  in  the  0=40"  laminate  is  higher  near  the  free  edge  then  in  the  0=15° 
laminate.  This  correlated  well  with  the  experimental  data  showing  higher  delamination 
initiation  loads  in  the  0=1 5"  plies. 

Multiscale  analysis  solution  and  calculated  bounds  for  the  first  and  second 
invariants  of  the  strain  tensor  in  the  matrix  phase  are  shown  on  Figure  7.  In  the  example 
at  hand  we  are  less  interested  in  fiber  characteristics,  since  the  initial  failure  mode  is 
matrix  cracking  and  delamination.  Figure  7  displays  the  values  of  the  first  and  second 
invariants  of  the  strain  tensor  calculated  at  the  midsurface  in  the  x=0  cross  section  as  a 
function  of  distance  from  the  hole  edge,  Figure  7. a  displays  results  obtained  by  using  ply 
level  analysis  and  Figure  7.b  results  obtained  by  using  SCSA  for  the  two  laminates.  In 
this  case  both  solutions  yield  qualitatively  similar  results.  Significantly  larger  values  of 
the  micro  level  stress  bounds  reflect  the  amplitude  of  the  strain  oscillation  in  the  unit  cell. 


00 

o  ° 
o 


8 


0  030 

0  025 

0  020 

0015 

0  DIO 

0  005 


0  000 


-0  52  *050  -0.48  *0  46  *0  44  *0  42  *0  40  -0  38 

V 

b) 

0.06 

0.05 

0  04 

0  03 


0  02 

0  01 

000 


-0.52  -0  50  -0  46  -0.46  -0  44  -0.42  -0.40 

Y 

Figure  7.  Comparison  of  the  distributions  of  strain  tensor  invariants  predicted  in  the 
matrix  phase  at  the  midsurface  at  x=0  as  a  function  of  distance  from  the  free  edge  by 
using  the  ply  level  solution  (a)  and  (b)  the  bounds  of  strain  tensor  invariants  predicted  by 
SCSA. 

2.2.  Normal  open-hole.  In  this  task,  we  compared  the  distribution  of  the  level  sets  of 
matrix  phase  strain  invariants  calculated  by  using  the  multiscale  solution  and  the  level 
sets  of  the  strain  invariants  calculated  from  the  macro  level  strains.  We  also  compared  the 
normalized  matrix  phase  strain  invariants  calculated  by  using  the  multiscale  solution  and 
the  level  sets  of  normalized  macro-level  stresses  and  their  combinations.  Such  a 
comparison  is  motivated  by  the  existing  approach  to  failure  prediction  in  composite 
laminates.  Thus,  the  failure  criteria  applied  on  the  micro  level  represent  the  Strain 
Invariant  Failure  Theory  (SIFT)  [6]  and  consist  of  the  first  and  second  invariants  of  the 
strain  tensor.  On  the  ply  level,  however,  the  most  common  failure  criteria  are  applied  in 
the  stress  space.  We  have  calibrated  both  failure  criteria  with  the  same  experimental 
values. 


Multiscale  solution 


a) 
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Comparison  of  invariants  of  the  homogenized  strain  field  to  those  of  the  macro  level 
strain  invariants.  The  first  invariant  bound  and  second  invariant  bound  will  be  denoted 
£M(e(z/’))  and  nM(e(u0)),  respectively.  A  0°  unidirectional,  4  in.  x  1  in.  rectangular  plate 
with  a  0.25  in.  diameter  open  circular  hole  was  considered.  Ply  thickness  was  0.005  in. 
A  1%  tension  strain  in  the  x-direction  was  applied  as  a  displacement  condition  ( ux  =  0.04) 
on  the  surface  x  —  2.  In  addition,  the  ^-displacement  was  restricted  on  the  surface  x  =  2. 
Similarly,  the  ^-displacement  and  ^-displacement  were  restricted  on  the  surface  x  =  -2. 
Midplane  symmetry  was  provided  for  the  laminate  by  the  boundary  condition  «.  =  0  on 
the  surface  2  =  0.  Thus  the  results  correspond  to  a  2-ply  model, 


figure  8(a)  and  8(b)  display  the  level  sets  of  the  second  invariant  of  the 
homogenized  strain  field  and  the  bounds  from  the  multiscale  solution  nM(e(t/0)), 
respectively,  figure  8  also  suggests  that  the  actual  strain  within  the  composite  can  be 
much  larger  than  the  state  predicted  by  the  homogenized  strain  field. 


(a)  Invariant  of  the  homogenized  strain  field 


(b)  Macroscopic  failure  criteria  nM(e(w0)) 


FiGURi-;  8:  Second  strain  invariant  comparison  for  unidirectional  0°  laminate 


Effect  of  thermal  prestress.  The  model  for  the  thermal  prestress  problem  was  identical 
to  the  model  of  the  previous  problem,  except  that  we  used  a  4-ply  laminate  (the  z- 
symmetry  condition  corresponds  to  an  8-ply  symmetric  laminate).  Each  ply  in  the 
laminate  had  a  thickness  of  0.005  in.,  and  thus  the  8-ply  specimen  was  0.04  in.  thick  in 
the  z-direction.  The  8-ply  laminate  was  an  IM7/977-3  [0/45/-45/90]s  composite.  The 
symmetry  plane  was  the  surface  z  =  0.  The  strain  invariant  bounds  in  the  matrix  phase 
were  computed  at  four  z-heights  in  the  model,  corresponding  to  the  middle  (with  respect 
to  z)  of  each  ply.  The  first  strain  invariant  in  the  45°  ply  is  shown  in  Figure  9.  The 
prestress  appears  to  increase  the  invariant  value  by  one  order  of  magnitude  on  nearly  the 
entire  ply.  This  trend  is  representative  of  the  prestress  effects  in  the  0°  and  ±45°  plies. 
The  prestress  did  increase  the  first  invariant  value  for  the  90°  ply,  but  by  only  a  few 
tenths. 

The  values  of  the  bounds  for  the  second  nM(e?(2z0))  in  the  45°  ply  are  displayed  in 
figure  10.  It  is  again  clear  from  the  figures  that  the  prestress  increases  the  second 
invariant  values  on  nearly  the  entire  plate.  The  prestress  does  not,  however,  appear  to 
significantly  increase  the  maximum  of  the  second  invariant  inside  the  ply.  It  is  seen  that 
the  prestress  increases  the  volume  of  the  sets  experiencing  the  higher  strain  states.  This 
volume  can  be  crucial  in  failure  predictions,  and  thus  its  increase  is  significant,  It  is 
important  to  note  that  second  invariant  values  in  the  ±45°  plies  are  about  twice  as  large  as 
the  values  in  the  0°  and  90°  plies.  The  off-axis  fiber  orientation  angles  clearly  experience 
larger  shear  strains. 
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ll  is  important  to  note  that  the  largest  strains  in  each  ply  are  concentrated  very 
near  the  hole  edge.  Hence  contours  approaching  these  extreme  values  are  not  visible  in 
these  figures.  The  invariant  contours  plotted  show  the  general  trends  of  the  invariants  in 
each  ply,  but  do  not  represent  the  extreme  invariant  values  that  occur.  For  the  first 
invariant,  the  largest  level-line  plotted  in  each  figure  is  roughly  50%  of  the  maximum 
invariant  in  the  ply.  For  the  second  invariant,  the  largest  level-line  shown  in  each  set  of 
figures  is  about  20%  of  the  maximum  invariant  in  the  ply. 


(a)  £M(e(w0))  without  prestress  (b)  £M(e(u0))  with  prestress 


FIGURE  9:  Macroscopic  failure  criteria  £M(e(«0))  in  the  +45°  ply 


(a)  nM(e(w0))  without  prestress  (b)  nM{e(w°))  with  prestress 

Figure  10:  Macroscopic  failure  criteria  nM(e(t/J))  in  the  +45°  ply 


Comparison  of  cubic  RVE  and  hexagonal  RVE.  We  created  a  hexagonal  unit  cell 
model.  Figure  5  shows  a  schematic  of  the  cubic  and  hexagonal  unit  cell  models, 
respectively.  The  macroscopic  computer  model  for  this  problem  was  an  8-ply,  symmetric 
laminate  with  an  open  hole,  and  was  identical  in  dimension  and  loading  to  the  model  for 
the  previous  problem.  Data  was  taken  with  and  without  prestress,  and  the  bounds  of  the 


(a)  Cubic  Unit  Cell  (b)  Hexagonal  Unit  Cel! 


Figure  !  1 :  Schematic  RVE  models  of  two  fiber  reinforced  geometries 

Figure  1 2  shows  the  first  invariant  bounds  £M(e(u0))  with  and  without  prestress. 
The  figure  suggests  once  again  that  the  prestress  increases  the  maximum  value  of  the 
macroscopic  failure  criteria,  as  well  as  increasing  its  value  on  the  entire  plate. 
Comparing  Figure  9  and  Figure  12,  it  is  also  evident  that  the  strain  state  is  smaller  when 
the  microstructure  is  hexagonally  periodic. 


(a)  £M(e(ua))  without  prestress 


■ 


(b)  £M(e(«0))  with  prestress 


Figure  12:  Bounds  for  the  first  strain  invariant  XM(e(u0))  in  the  +45°  ply 


We  note  here  that  the  color  scales  are  the  same  for  Figure  9  and  Figure  1 2,  as  well 
as  for  Figure  10  and  Figure  13.  Thus  the  corresponding  figures  for  the  cubic  RVE  and 
the  hexagonal  RVE  can  be  compared  directly  by  observing  the  color.  Figure  13  plots 
second  strain  invariant  bound  UM(e(«0))  with  and  without  prestress.  Once  again,  the 
macroscopic  failure  criteria  values  calculated  using  the  hexagonal  RVE  were  less  than  the 
values  calculated  using  the  cubic  RVE. 


(a)  nM(e(i(0))  without  prestress  (b)  nM(e(«0))  with  prestress 


Figure  13:  The  bound  for  the  second  strain  invariant  nM(e(M0))  in  the  +45°  ply 

Failure  analysis  of  a  symmetric  laminate  with  an  open  hole.  We  compared  the 
normalized  bounds  of  the  matrix  phase  strain  invariants  with  the  normalized  macro-level 
stress  values  and/or  their  combinations.  Experimental  strength  data  were  used  to  back- 
calculate  critical  strain  invariant  values  for  the  matrix  phase  of  the  composite.  These 
critical  values  were  used  to  normalize  the  strain  invariant  bounds.  Hashin  criteria  and 
maximum  stress  criteria  were  applied  to  the  ply  level  (homogenized)  stresses.  The  over¬ 
stressed  (over-strained)  zone  is  the  volume  of  the  composite  in  which  the  value  of  the 
failure  criteria  meets  or  exceeds  any  of  the  critical  strength  parameters  (indicated  by 
yellow  in  Figures  1 5  -  19). 

The  computer  model  for  this  experiment  was  identical  in  dimension  to  the  model 
of  the  previous  two  experiments.  The  loading,  however,  was  reduced.  Since  complete 
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loss  of  load  carrying  capacity  typically  occurs  between  0.6%  strain  and  0.7%  strain,  wc 
computed  over-stressed  zones  at  0.3%  and  0.6%  uniaxial  tension  (applied  via 
displacement  on  the  surface  „y  =  2).  All  other  boundary  conditions  were  the  same  as  in 
the  prior  experiments. 

figures  14.  15.  and  16  display  the  overstressed  zones  for  the  phase  strain  invariant. 
I lashin  failure  criteria,  and  maximum  stress  failure  criteria,  respectively  at  0.3%  tension 
strain  loading.  At  such  a  small  loading,  the  0°  ply  has  very  small  overstressed  zones  just 
at  the  hole  surface.  The  45°  ply  shows  moderately  small  over-stressed  zones  protruding 
from  the  hole  surface  out  toward  the  ply  edge.  The  -45°  ply  (not  shown)  was  similar  to 
the  45°  ply.  and  the  90°  ply  (also  not  shown)  was  almost  completely  overstressed  at  just 
0.3%  strain.  The  three  failure  theories  agree  quite  well  on  the  location  and  size  of  the 
over-stressed  zones.  Moreover,  the  theories  appear  to  agree  on  the  failure  mode.  The 
dashed-line  boundary  around  the  macroscopic  failure  criteria  overstrained  zone  (Figure 
14(a))  indicates  shear  failure.  The  solid-line  boundary  around  the  I  lashin  overstressed 
zone  (Figure  13(a))  indicates  failure  via  a  quadratic  combination  of  transverse  tension 
and  in-plane  shear  forces.  The  boundary  around  the  maximum  stress  over-stressed  zone 
(Figure  14(a))  indicates  failure  by  the  in-plane  shear  force.  Thus  the  results  suggest  that 
the  three  theories  agree  that  the  shear  forces  will  initiate  failure  at  the  hole  edge  in  the  0° 


ply, 


Figures  17,  18.  and  19  visualize  the  over-stressed  zones  for  the  micro  level  strain 
bounds,  Flashin  failure  criteria,  and  max  stress  failure  criteria,  respectively  at  0.6% 
tension  strain  loading.  The  90°  and  ±45°  plies  all  have  large  over-stressed  zones  (only 
+45°  ply  is  shown).  The  small  overstressed  zones  very  near  the  hole  edge  in  the  0°  ply- 
are  the  shear  overstressed  zones  seen  at  0.3%  strain.  These  shear  stresses  near  the  hole 
edge  produce  the  characteristic  cracks  running  parallel  to  the  fibers  in  the  0°  ply.  The 
slightly  larger  overstressed  zones  in  the  0°  ply  appear  to  be  transverse  tension  failure. 
These  overstressed  zones  are  believed  to  be  an  effect  of  the  ply  stacking  sequence.  Once 
again,  the  micro-level  strain  invariants  based  failure  criteria  appear  to  be  consistent  with 
the  predictions  of  the  Hashin  failure  criteria  and  the  maximum  stress  failure  criteria.  In 
addition  to  the  three  different  failure  criteria  being  in  agreement  on  the  size  and  location 
of  the  overstressed  zones,  they  also  once  again  appear  to  agree  on  the  failure  mode  in 
each  case. 
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(b)  +45°  Ply 


FIGURE  1 8:  Hash  ill  failure  criteria  overstressed  zone  at  0.6%  strain 


(a)  0°  Ply  (b)  +45°  Ply 


Figure  19:  Maximum  stress  failure  criteria  overstressed  zone  at  0.6%  strain 


TASK  3.  In  this  task  we  tried  to  understand  the  failure  sequence  of  the  composite  scarf 
repair  patch  and  construct  a  meaningful  goal  function  for  strength  retention  optimization. 
The  main  development  in  this  task  was  to  model  nonlinear  deformation  and  failure  in  the 
adhesive.  The  deformation  of  the  composite  and  the  patch  was  modeled  with  less  detail 
on  ply  level.  We  heavily  utilized  experimental  data  obtained  by  Capt.  B.M.  Cook  [8], 
which  included  fabrication  and  testing  of  the  scarf  repair  as  shown  on  Figure  20. 


FIGURE  20.  Experimental  setup  for  tensile  testing  (a)  and  a  plane  view  of  the  failed  scarf 
repaired  specimen  (b). 

The  key  Finding  in  this  task  was  that  modeling  of  the  failure  of  the  adhesive  in  the 
nonlinear  regime  is  crucial  to  predicting  the  strength  of  the  repair.  We  found  that  despite 
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the  low  stiffness  the  adhesive  transfers  enough  load  in  the  elastic  regime  to  practically 
remove  the  stress  concentration  due  to  scarfing  and  patching  thus  leading  to  95%  strength 
retention  prediction.  It  is  the  adhesive  failure  that  drastically  increases  the  stress 
concentration  in  the  parent  panel  and  leads  to  premature  failure.  Spline  approximation 
based  3D  ply  level  analysis  was  used  for  failure  prediction.  The  plate,  adhesive  and  the 
repair  patch  models  are  shown  on  Figure  2 1 . 

Patch 


Figure  2 1 .  Computational  model  for  the  scarfed  plate,  patch  and  adhesive. 

A  nonlinear  elastic  model  was  utilized  to  model  the  adhesive  At  the  same  time  the  fiber 
failure  in  the  scarfed  panel  was  predicted  by  using  the  critical  failure  volume  (CFV) 
method  recently  proposed  in  [9-10],  The  average  applied  load  for  fiber  failure  was 
predicted  at  each  load  level  characterized  by  a  state  of  delamination  in  the 
adhesive/scarfed  panel  interface.  Initially,  the  stress  distribution  in  the  repaired  panel  is 
uniform  and  the  fiber  failure  stress  predicted  is  identical  to  that  of  virgin  laminate.  The 
increase  of  the  applied  load  leads  to  initiation  and  growth  of  the  delamination  between 
the  repair  patch  and  the  plate.  The  stress  distribution  in  the  patch  becomes  nonuniform 
with  the  stress  concentration  typical  of  open  hole  composites.  At  some  applied  load  the 
average  predicted  fiber  failure  load  and  the  applied  load  become  equal.  This  load  value  is 
considered  the  failure  load  for  the  repaired  panel. 

The  unidirectional  properties  of  the  1M6/3501-6  material  were  as  follows: 
E|  i  =  l  75.27Gpa,  E22=E33=9,79Gpa.  v,3=  v13=0.33,  v23=0.49,  G)2=Gi3=5.5  1  Gpa  and 
G23=2.96Gpa,  with  the  coefficient  of  thermal  expansion  of  an=4.10'6C°  and  a22-32  10' 
6C°  The  Weibull  parameters  for  the  strength  in  the  fiber  direction  were  equal  to 
Xt=2.06Gpa,  V0=  1638.7mm3  and  a=40.  The  value  of  the  Weibull  modulus  was  taken 
from  Wisnom  et  al.  [1 1  ].  The  /nim=0.266mm  was  calculated  according  to  reference  [9-10]. 
The  initial  elastic  properties  of  the  adhesive  were  E=3.1Gpa,  v=0.38  and  the  thermal 
expansion  coefficient  was  a=0.7  I0'6C°.  The  ply  and  adhesive  thickness  were  h=0.!3mm. 
The  nonlinear  shear  deformation  curve  t(y)  for  the  FM-300M  adhesive  was  provided  by 
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Structural  Integrity  Branch  of  the  Air  Force  Research  Laboratory’s  Materials  and 
Manufacturing  Directorate  [12]  and  is  shown  in  Figure  21.  Three-  dimensional  nonlinear 
deformation  of  the  adhesive  was  describer  by  assuming  nonlinear  elasticity  with  constant 

poisson  ratio  and  shear  modulus  equal  to  G=r(2/V3  svm).  where  eVm-  is  Von-Mises  strain, 
where  G(y)=dx(y)/dy.  The  temperature  change  AT=- 1 55°C. 


Figure  2 1 .  Nonlinear  shear  deformation  curve  of  the  for  the  FM-300M  adhesive. 


At  low  loads  prior  to  the  adhesive  failure  the  stress  field  in  the  scarfed  panel  is  quite 
uniform,  as  shown  on  the  Figure  22.  No  stress  concentration  typical  of  an  open  hole  is 
present.  The  average  fiber  failure  load  is  predicted  equal  to  that  of  virgin  panels. 

(a)  b) 


Figure  22.  Normal  oxx  stress  distribution  (a)  in  the  adhesive  and  (b)  the  scarfed  plate  prior 
to  any  dclamination  on  the  adhesive/panel  interface 
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In  the  case  of  scarfed  open  hole  panels  a  severe  stress  concentration  typical  of  an  open 
hole  exists,  Figure  23.a.  The  strength  of  such  panels  is  about  1/3  of  the  virgin  panels  and 
is  predicted  accurately  by  using  the  CFV  method  [9-10].  Thus  the  strength  of  the  repair  is 
determined  by  the  adhesive  failure  and/or  its  nonlinear  behavior.  The  stress  distribution 
in  the  parent  panel  corresponding  to  partially  failed  adhesive  is  shown  on  Figure  23b. 
This  stress  concentration  is  less  severe  than  in  the  case  of  open  scarfed  hole  and  the 
predicted  strength  values  of  the  repaired  panels  are  at  the  level  of  2/3  of  the  initial  virgin 
panel  strength. 


Figure  23.  Normal  ctxx  stress  distribution  in  the  scarfed  and  repaired  panels,  (a)  Open 
scarfed  hole  and  (b)  partially  delaminated  adhesive  at  90%  of  the  scarf  failure  load. 


Experimental  data  and  predicted  strength  values  for  small  and  large  size  unnotched 
tensile  coupons,  as  well  as  strength  of  scarfed  open  hole  and  scarfed  repaired  panels  are 
shown  on  Figure  24.  Significant  reduction  in  experimentally  measured  strength  of  large 
virgin  panels  as  compared  to  small  coupons  is  attributed  to  grip  failure  of  most  large 
virgin  specimens.  Only  5%  strength  reduction  due  to  volume  increase  is  predicted 
theoretically.  All  unrepaired  and  repaired  scarfed  laminates  exhibited  brittle  failure 
through  the  cavity,  thus  rendering  the  test  data  valid  for  comparison  with  theoretical 
predictions.  The  predicted  strength  values  for  both  cases  are  close  to  those  measured 


Figure  24.  Experimental  and  numerical  results 
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experimentally. 

TASK  4.  Development  of  an  automated  mutiscalc  mechanism  within  the  B-splinc 
Approximation  Method  (BSAM)  computer  program  was  accomplished  ahead  of  schedule 
and  utilized  to  produce  all  results  demonstrated  in  this  report. 

TASK  5. 

In  the  course  of  the  project,  we  have  understood  that  our  original  vision  of  repair 
optimization  based  on  elastic  properties  of  adhesive  was  too  simplistic.  Indeed  as  seen 
from  the  previous  task  the  repair  strength  is  a  direct  function  of  the  adhesive  failure, 
w'hich  determines  the  load  transfer  between  the  repair  patch  and  the  parent  laminate.  An 
optimization  algorithm  which  would  be  based  on  a  repair  strength  as  a  goal  function  will 
be  highly  impractical,  requiring  nonlinear  analysis  at  every  optimization  iteration  and 
preventing  the  use  of  gradient  methods  [13]  for  optimization.  These  complications 
delayed  our  progress  and  required  to  propose  alternative  optimization  ideas.  The  one  we 
have  selected  for  implementation  is  to  minimize  an  elastic  Lpnorm  of  a  stress  invariant 
failure  criteria  in  the  adhesive.  Such  an  approach  is  targeted  to  the  delay  of  failure 
initiation  in  the  adhesive.  Although  not  completed,  we  have  obtained  encouraging  results. 
Recall  figure  24,  where  the  predicted  strength  of  the  repaired  composite  was  at  60Ksi. 
This  result  was  obtained  for  a  traditional  ply  by  ply  replacement  composition  of  the  repair 
patch.  After  we  switched  two  ply  orientations,  namely  0(l  and  90,)  inside  this  patch,  we 
delayed  the  onset  of  adhesive  failure  to  65Ksi.  We  expect  to  further  delay  the  adhesive 
failure  by  application  of  the  optimization  algorithm.  These  developments  are  being 
closely  followed  by  the  Structural  Integrity  Branch  of  the  Air  Force  Research 
Laboratory’s  Materials  and  Manufacturing  Directorate,  which  is  presently  conducting  an 
experimental  program  to  verify  our  findings. 

CONCLUSIONS 

1 )  We  integrated  the  rigorous  microlevel  (fiber,  matrix)  stress  bound  estimates  by 
using  simply  coupled  scale  analysis,  with  the  spline  approximation  stress  analysis  tools  * 
BSAM.  A  robust  multiscale  analysis  framework  was  developed  and  applied  within  and 
beyond  the  scope  of  the  present  grant. 

2)  A  key  development  was  the  extension  of  the  stress  bound  estimates  in 
heterogeneous  materials  with  periodic  microstructure  to  include  regions  of  nonperiodic 
microstructure.  Such  regions  included  ply  interfaces  in  fiber  reinforced  composite 
materials,  which  are  present  in  all  laminated  composite  structures.  The  improved 
accuracy  of  the  multiscale  analysis  was  verified  by  comparison  with  three-dimensional 
analysis  of  model  cross  ply  composite  laminate,  where  fibers  were  modeled  explicitly. 

3)  The  developed  capabilities  were  applied  to  composite  scarf  repair  optimization. 
Our  analysis  addressed,  for  the  first  time,  the  composition  of  the  repair  patch  in  addition 
to  more  conventional  geometric  parameters,  such  as  scarf  angle  and  repair  depth.  We 
found  that  altering  the  conventional  ply  by  ply  replacement  schema  of  the  repair  patch 
does  not  lead  to  premature  failure,  but  delays  the  failure  of  the  adhesive  and  increases  the 
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strength  of  the  repaired  composite.  This  work  is  coordinated  with  Structural  Integrity 
Branch  of  the  Air  Force  Research  Laboratory’s  Materials  and  Manufacturing  Directorate, 
which  is  presently  executing  an  experimental  program  to  verify  our  findings. 

4)  We  have  transitioned  the  periodic  RVE  based  micromechanical  failure  criterion 
evaluation  algorithm  to  Sikorsky  Aircraft,  where  it  was  applied  to  failure  analysis  of  a 
helicopter  flex  beam.  We  are  also  applying  our  methodology  to  investigate  applicability 
of  various  rnicromechanical  failure  criteria  to  failure  prediction  of  composites  under 
complex  states  of  stress. 
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